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The fluid forces on a centrifugal impeller rotating and whirling in a vaned 
diffuser are analysed on the assumption that the number of impeller and diffuser 
vanes is so large that the flows are perfectly guided by the vanes. The flow is 
taken to be two dimensional, inviscid, and incompressible, but the effects of 
impeller and diffuser losses are taken into account. 

It is shown that the interaction with the vaned diffuser may cause 
destabilizing fluid forces for whirling motion even at design flow rate. 
Discussions are made to elucidate the physical mechanisms that produce the 
destabilizing fluid forces. From these discussions, it is found that the whirling 
forces are closely related to the steady head-capacity characteristics of the 
impeller. This physical understanding of the whirling forces can be applied also to 
the cases with volute casings. 

At partial capacities, it is shown that the impeller forces change greatly when 
the flow rate and whirl velocity are near to the impeller or vaned diffuser 
attributed rotating stall onset capacity, and the stall propagation velocity, 
respectively. In such cases the impeller forces may become destabilizing for 
impeller whirl. 


INTRODUCTION 

It has been made clear by pioneering experiments by Ohashi et al. (ref. 1) and 
Acosta et al. (refs. 2 and 3) that the fluid forces on centrifugal impellers may 
become destabilizing for whirling motion under certain conditions. Even if the 
impeller forces are stabilizing, it is important to estimate the forces since there 
are many cases in which the whirling instability is caused by bearing or seal forces. 

So far the following examples have been reported in which the impeller forces 
become destabilizing for whirl. Ohashi et al. (ref. 1) have measured whirling 
forces on two-dimensional impellers in vaned or vaneless diffusers to capture the 
basic nature of the forces and to compare with their two-dimensional potential flow 
analyses (refs. A and 5). In this, it was found that the forces are basically 
stabilizing but become destabilizing for whirling with small positive whirl velocity 
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ratios at partial capacity. Recently new testing equipment was developed (ref. 6), 
and preliminary results on a three-dimensional boiler feed pump impeller model 
whirling in a vaned diffuser were reported. The impeller forces become very large 
compared with the 2D cases, but they are stabilizing both at design and shut off 
capacities . 

Acosta et al. focused on the interaction with the volute casing and have shown 
that the forces on impellers whirling in a volute may become destabilizing even at 
design capacity (refs. 2, 3 and 7). A two-dimensional unsteady flow analysis 
(ref. 8) was carried out assuming that the impeller has an infinite number of vanes, 
and it was found that it can simulate the experimental results fairly well. 

Bolleter and Wyss (ref. 9) used a "rocking arm" test apparatus to measure force 
matrices of a boiler feed pump impeller with outlet guide vanes. Whirling forces 
were synthesized on the assumption of the linearity of the disturbance, and it was 
shown that the forces are destabilizing for small positive whirl even at design 
capacity. 

Two of the present authors have made an analysis (ref. 10) of impeller forces 
on a two-dimensional impeller whirling in a vaneless diffuser. There, it was shown 
that the whirling forces are basically stabilizing, as shown by Ohashi et al. 
(ref. 1), but they may become destabilizing if the flow rate and the whirling 
velocity are near the onset flow rate and the stall propagation velocity of rotating 
stalls due to the impeller or the vaneless diffuser. 

The present study treats the cases when a impeller whirls in a vaned diffuser 
on the assumption of a two-dimensional inviscid incompressible flow. The number of 
impeller and/or diffuser vanes is taken to be infinite, and the total pressure 
losses in the impeller and the diffuser are taken into account. It is assumed that 
the outlet of the diffuser is open to a space of constant pressure. The problem is 
linearized on the assumption that the eccentricity and hence the unsteady components 
are sufficiently small. 


FLOW ANALYSIS 

We consider the case when, as shown in Fig. 1, a logarithmic centrifugal 
impeller with a vane angle (5 and inner and outer radii r^ and r£ rotates with 
an angular velocity Q and whirls with an eccentricity e and an angular velocity 
a), in a logarithmic vaned diffuser with a vane angle a, inner and outer radii r 3 
and r 4 * We take a stationary frame z' * x' + iy' = r'e 1 ®' with its origin at the 
center of the diffuser, and a moving frame z = x + iy = re*® with its origin at 
the center of the impeller translating with its axes parallel to the stationary 
frame. We put ' on the velocities relative to the z* frame, and the velocities 

without 1 represent the velocities relative to the z frame. The total velocity 

and the steady and unsteady components are represented by v, V, and v^, 

respectively. For unsteady components, we use a complex representation with respect 
to an imaginary unit j such as ^rpd = ^rp exp{j(wt - 9)} where the suffix r 

means radial component and p = 1, 2, 3, k means the values at r = r^, r 2 * r 3 , and 
r 4 * We consider that the real part carries the physical meaning. 

Relations governing unsteady components 

The results of the relations and expressions obtained for the vaneless diffuser 
case (ref. 10) could also be used here. 
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Velocity at the impeller inlet: The velocity disturbance is composed of that 
due to relative dislocation of a source Q at the origin of the z* frame and an 
unknown potential disturbance 

v r i = -eQ/(2irri 2 ) + v rl * (1) 

VQi = -eQj/(2irr 1 2 ) - jv rl * (2) 

where v r i* is an unknown complex constant. 


Velocities at the impeller outlet and the diffuser inlet (p = 2, 3): The 

steady velocity in the vaneless part (r 2 < r < r^) is given by (V r , Vq) = 
(Q/2irr, I72irr) where I = “ Q/ ( 2TTr 2 ) cot(3] is the impeller circulation. 

We assume a vorticity C shed from the impeller is transported at this steady 
velocity. Then the unsteady velocities can be given as a sum of the velocity 
induced by the vorticity field of magnitude ^2 anc * a potential flow disturbance of 
strength A. 


where 


<1 

II 

Z r pC2 + A rp A 




(3) 

V0P = 

Z© p C2 + A0p A 




(4) 

2 rp = 

(Rr " jRl>rp 

- ( r 2/ r p) 2 e 

-2j3 

(Rr - jRl ) r 2 


2 0p = 

(®R * J®l)rp 

- j (r 2 /r p > 2 

e- 2 jP 

(Rr - jRj) r 2 


A rp = 

1 - (r 2 /r p ) 2 

e -2 ^, Ae p = 

-jU 

+ (r2/r p ) 2 e -2 ^} 



Rr - jRl =-j(S + T), ©r- j®i = S - T 


S(r ) = (1/2) [ Mrg, 0)(rQ/r) 2 drg 
T(r) = (1/2) J 3 X(r 0 , 0) dr 0 


\(r,0) = exp( [-j{irto(r 2 - r 2 2 )/Q + 9 - (r/Q)log(r/r2>} ] 

The expressions (3) and (4) are so determined that the flow tangency condition at 
the impeller outlet v @2 = -v r 2 cot|3 is satisfied. 

Continuity of mass flow between impeller inlet and outlet: Defining 4> = ©]_ — 

02 = cot3 log(r 2 /r^) as shown in Fig. 1, we obtain the following continuity relation 
between r = rp and r 2 * 

r iv ri = r 2 v r 2ej ( t> (5) 

Strength of shed vorticity: We obtain the following expression by integrating 

the momentum equation from inlet to outlet of the impeller and using the 0 
component of Euler's equation at the impeller outlet. 
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r 2 v r2 ^2 = J ^ r 2 + D 


( 6 ) 


where 

J = (to - Q){r 2 (l + j cot|3) + M} -jLi 
D = cQe~J < t ) {2r^(to - Q) - L 2 } /(2irr].2) 

' r 2 _ 9 

and M = r 2 /(r sin 3 )ds = ^(sinft) dogt^/r^) is the effective length of the flow 

r 3 

channel of the impeller, and L 2 are parameters representing the effects of 

total pressure loss in the impeller. These are given as follows. We assume that 

the loss can be given as a sum of incidence and hydraulic losses: 

Ap t /p = U s (l/t a - l/t ad )2 + £*}v rl 2/2 = Z(t a )v rl 2/2 (7) 

where t a is the tangent of the mean flow inlet angle (3 ^ » and t ad = tan|3. The 

hydraulic and incidence loss coefficients and £ s are evaluated from 

C* = 2(1 - rih)( r l/ r 2) 2 ( 1 “ *f*d cot(3 )/$ d 2 

C s = 0.3 + 0.63°/90° 

using the shock free flow coefficient $ d and the impeller efficiency % at $ d . 
For simplicity, we neglect the delay of loss and evaluate the loss from quasi-steady 
application of equation (7). After linearization on the assumption that the 

disturbance is small, we obtain the following expressions for and L 2 : 

L 1 = V rl (r 2 /r 1 ){Z + (l/2)tan3i(l - j tan^Z'} 

L 2 = V rl Z ' tan2(3i 

The condition at diffuser outlet 

We assume that the outlet of the diffuser is connected to a space of constant 
pressure. From this condition, we derive a relation that should be satisfied at the 
diffuser inlet. Taking the total pressure loss Ap td in the diffuser into account, 
we can express the unsteady Bernoulli's equation as follows. 

f r 4 

P 3 /p = P 4 /P + 1 / 2 (v 4 ' 2 - V 3 * 2) + 3v s ’/3t • ds + Ap td /p ( 8 ) 

r 3 

We evaluate Ap td by equation (7) after replacements of Z(t a ) with Z d (tan a 3 ) 
and v r i with v’ r 3 . After the linearization of equation ( 8 ) with respect to the 
unsteady component and the differentiation with 0 ', we obtain 

1/p • 3 p 3 / 3 © 3 ' = -j[(E d + Li d )v r 3 ' - (I72irr3 + L 2 d / 2 )vQ 3 '] (9) 

where 

E d = V r 4 (r 3 /r 4 )(l + cot 2 a ) - V r 3 + jtoM d 
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M d = effective length of diffuser flow channel 
= r3(sina) -2 log(r3/r4) 

Lid = ZdVr 3 + (l/ 2 )V r3 Z d ’ tani 3 
L 2 d = v r 3 z d' tan2(X 3 

and the relation 9p4/903’ = 0 has been used. 013 = tan - l(Q/r) is the mean flow 
angle at the diffuser inlet. 

The unsteady component of the absolute velocity at the diffuser inlet is given 
as a sum of the unsteady component in the moving frame, a component due to the 
dislocation of the steady flow field (V r , V©), and the whirling velocity: 

v r 3' = v r 3 + e(Q - jr)/2irr3 2 ) + ewj 

v@3 ' = VQ3 + e(T + jQ)/2irr3 2 ) + eto 

Putting equations ( 9 ) and ( 10 ) in the 0 ' component of a linearized Euler’s 
equation, we finally obtain the following equation representing the diffuser outlet 
condition. 


P; C 2 + P A A = G 


(ID 


where 


P C = (j« + V r 3/r 3 + jL 2 d/ 2 r 3 )Z 0 3 - j(E d + L ld )Z r3 /r3 + V r 3Z a 

P A * (j« + V r 3/r 3 + jL 2 d / 2 r 3 )A 03 - j(E d + L ld )A r3 /r 3 - jZ b V r3 

G = -e(jw + V r 3/r3 + jL 2d /2r3)(w + V03/r3 + jV r 3/r3) 

+ ej(E d + L ld )(V r3 /r3 + jto - jV03/r 3 )/r 3 + 2eV r3 (V0 3 + jV r3 )/r 3 2 

Z a = X(r3, 0 ) - 2 S(r 3 )/r 3 + 2 j (r 2 2 /r 33 )e~ 2 J& (Rr - jRi) r 2 

Z b = - 2 (r 2 2 /r 33 )e -2 j^ 

Determination of flow and impeller forces 

From equations ( 1 ) to (6) and ( 11 ), we obtain 

A ■ {G(r 2 V r2 - JZ r2 ) - P^ D}/{P^ JA r2 + P A (r 2 v r 2 ~ JZ t2 )) ( 12 ) 

Other unknowns are determined by the successive use of equations ( 1 ) to (6) and 
( 11 ). Each term of the numerator of the above equation has a factor e. By 
considering the case with e =0, we find that the rotating stall onset conditions 
under the interaction of impeller and diffuser are obtained by setting the 
denominator of equation (12) to be zero. 

The impeller forces are determined from the balance of momentum. If we 
separate the impeller forces into components normal to (n) and tangential to (t; 
positive counter clockwise) the whirling orbit, each component stays constant. The 
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impeller forces (F n , F t ) normalized with pir(r 2 Q) 2 e are given by summing up the 
following components. (0 E Q/2Trr 2 2 Q is a flow coefficient.) 

Force due to inlet pressure 

(F n + jF t )pi = (ri/r2)[(a)r^/Qr2){-2j(r 2 /r 1 ) 2 (j> + co/Q} 

- v r 2 /(eQ)ej^{ jco/Q + 4> ( 3r 2 / ^1 ) 2 } ] (13.1) 


Force due to inlet momentum 

(F n + jF t ) m ]_ = 0 [ 2<J> ( r 2 / ) 2 + 2 jco/Q + v r2 (r 2 /r i)/eQ] (13.2) 

Force due to outlet pressure 

(F n + jF t ) p 2 = -(c o/Q) 2 - je'3^ (L 2 /r 2 ft)0(r 2 /ri) 2 

- 2 j 4 >(r 2 /r^ )e“3^(l - co/Q) 

- {j(l + M/r 2 )(l - co/Q) - 0(1 + cot 2 3) 

- L 1 /(r 2 (2)}v r2 /(eQ) 

Force due to outlet momentum 

(F n + jF t ) m2 = -{20 + j(l - 20 cot3}v r2 /(eQ) - 2j0(co/Q) 

Force due to change of momentum of the fluid in the impeller 
(F n + j F t)m = -j(u/Q){l - (r 1 /r 2 )ej<|)}0 r2 /(eQ) + (co/Q) 2{1 - ( ri /r 2 )2} (13.5) 

RESULTS AND DISCUSSIONS 

We consider an impeller with ~ 0-5* 3 = 30°. The shock free flow 

coefficient for this impeller is 0^ E v r2 /(r 2 Q) = 0.1443, which is looked upon as 
the "design" condition and is also shockless for diffusers with a = 11°. The loss 

coefficients of the impeller were estimated from the equations shown in the last 

section assuming Hh = 0*9. The hydraulic loss coefficient for the diffuser was 

estimated from the expression of the coefficient after multiplying by (r 2 /r^) 2 and 
using nh = 0*9. The diffuser incidence loss coefficient was obtained from the 
expression for the impeller after replacing 3 with a. 

Impeller forces at design flow rate 

Figure 2 shows the plots of F n and F t against the whirl velocity ratio 
co/Q. There, the diffuser inlet radius r 3 is changed keeping (r^ - r 3)/ r 2 t0 
0.25. It has been known that if F t is in the same direction as the whirl it is 
destabilizing, and vice versa. We find a destabilizing region for small positive 

whirl for r 3 / r *2 < 1*5. This region diminishes with the increase of r 3/ r 2* A 

similar destabilizing region is observed also for the cases with a volute. 

Discussions on this will be made shortly. F n for w/Q < 0 increases with the 
decrease in r 3 /r 2 * which is considered to be caused by the increase of the apparent 
mass due to interaction with the vaned diffuser. If we increase r 3/ r 2> F n anc * 
F t tend to the values for the case with a vaneless diffuser with outer radius 
r 3/ r 2 = 5.25. 


(13.3) 


(13.4) 
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Figure 3 shows the effects of r 4 /r£. We see that F n for o>/D < 0 increases 
with the increase in but is not affected by r^/^ systematically. 

Note that the destabilizing region exists even for the case with r^ = r 3 . On the 
other hand, if we change the vane angle a and keeping the effective length 

of the diffuser to be constant, as shown in Figure F n does not change much 

but F t for small positive to negative co/Q decreases with the increase in a. Now 

we can say that F n for co/Q < 0 is mainly affected by M^. F t for w/Q £ 0.4, 

and hence the existence of the destabilizing region, is largely dependent on the 
value of a. 

In Figure 5 the impeller forces are separated into components corresponding to 
equations (13.1) to (13.5). We see that the component due to outlet pressure is 
dominant for co/fi < 0.4. Although not shown here, the contribution of outlet 

pressure is larger for the cases with smaller ^*3 / r *2 • This is true also at partial 

capacities . 

In all of the above calculations the effects of impeller and diffuser losses 

are taken into account. The above tendencies are not changed by neglecting the 

losses. After all, it is considered that the flow turning effect of the vaned 

diffuser is responsible for the existence of the destabilizing region. 

Let us consider the physical mechanism with which the impeller forces become 
destabilizing. To simplify the problem we take the limits r 3 ■* r 2 and a 0. 

For such cases we can obtain a purely analytical solution without numerical 

integrations with respect to the shed vorticity. By letting a *♦ 0 in the 

continuity equation (v^ 1 )^ = (r 3 /r 4 )^(v r 3 ' )^/sin^a, we find that a small increase 
in v r 3 * some part of the circumference results in a large increase in V 4 * and 

hence in P 3 there, from equation ( 8 ). Then the outflow from the impeller to that 
part will be diminished; a diffuser with a small a has a smoothing effect on the 
circumferential flow rate distribution through the diffuser. By letting v r 3 * ■* 0 
and r 3 r 2 in equation ( 10 ), we obtain 

v r 2 * -e(Q - jD/(2irr2 2 ) - euj 

= eQ[-<j> + {(1 - <|> cotf3) - co/Q}j] (14) 

This means that an outflow distribution v r 2 from the impeller is so determined as 
to cancel the radial absolute velocity disturbances due to the dislocation of the 
flow field of impeller circulation I and source Q and due to the whirling 

velocity. This equation can be obtained also by letting a + 0 in the analytical 

solution. 

We estimate the impeller forces due to a pressure distribution which is 
obtained by assuming that the inertia on the fluid in the impeller due to v r 2 is 
balanced by the outlet pressure distribution. Then we have 

(F n + jF t ) p 2 = -j(l - w/Q)Mv r 2 /(r 2 eQ) (15) 

From the definition of complex representation, we find that the complex plot 
of v r 2 gives the azimuthal direction where the radial velocity becomes maximum and 
the amplitude of the radial velocity fluctuation, in a plane with its real axis in 
the direction of eccentricity (n) and the imaginary axis in the tangential (t) 
direction. Equation (14) shows that the radial velocity will become maximum in the 
direction (n,t) = (-$, ^ E 1 - $ cot|3) for w/Q = 0 and the pressure in the direction 
(-«|>, -<j>) where <J> and Y are flow and head coefficients. The resulting force is. 
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from equation (15), in the direction of (<{», 4>), i.e., outward and in the direction 
of rotation. Thus the impeller force is closely related to the steady impeller 
characteristics. F n and F t determined from equations (14) and (15) are given in 
the APPENDIX. Next is considered the effects of impeller loss. Near design flow 
rate, the loss becomes maximum in the direction of v r 2 - There the outlet pressure 
is diminished and the resulting force can be expressed as Z(r 2 /ri)^ <|> v r 2 /(eQ), 
which is in the same direction as v r 2 . 

Figure 6 shows the above relations for a •* 0 and *31*2 * 1 which impeller 
and diffuser losses are neglected. This figure shows that equation (14) can 
estimate v r 2 fairly well for a = 1° and that the impeller force may be 
approximated by the outlet pressure force of equation (15) for u/Q < 0. Figure 7 
shows an example for a more realistic case with *11*2 = 1*05, a = 11° with the 
effects of losses taken into account. The absolute value of v r 2 is significantly 
smaller than that estimated from equation (14), which is mainly caused by the 
finiteness of a. We see that the impeller force may be estimated from the simple 
relations if we use the exact value of v r 2 in equation (15) and take into account 
the effect of loss discussed above. The real part of v r 2 , and correspondingly F^ 
for w/Q = 0, approaches zero by increasing t^/t2 and hence the destabilizing 
region diminishes. 

From the above discussions we may understand the mechanisms of the 
destabilizing impeller forces as follows. A vaned diffuser with a small vane angle 
a has the effect to smooth out the circumferential flow rate distribution through 
the diffuser. This effect brings about a circumferential flow distribution through 
the impeller to cancel out the absolute radial velocity disturbances at the diffuser 
inlet due to the dislocation of the flow field of impeller circulation I and 
source Q, and to the whirling velocity. Such a flow distribution through the 
impeller is established by a pressure distribution around the impeller, which 
balances with the inertia on the fluid in the impeller caused by the flow 
distribution. In cases of a •* 0 and ^ 3/^2 * 1 the flow through the impeller for 
co/Q = 0 becomes maximum in the direction of (n,t) = (~4>,<J/) and the impeller force 
(F n , F^) due to the outlet pressure distribution caused by the flow distribution is 
positively proportional to (q/ , <J>). This tendency is carried over to general cases 
of finite but small values of a and *"$1*2 * 1, resulting in a destabilizing 
region at small positive u/Q. 

The tendencies that the outlet pressure force is more dominant for the cases of 
smaller *z/*2 an< * ®/Q ma y be understood by considering that the above mechanism 
becomes more definitive for those cases. When the diffuser angle is not constant, 
the outlet angle will have a more important effect on the impeller forces. In cases 
with a volute, the radial velocity on the volute is strictly controlled by the 
volute surface. Hence the latter half of the above discussions can be applied also 
for the cases with a volute. 

Impeller forces at partial capacity 
(Cases with infinitesimally small impeller/diffuser gap) 

We consider the cases with ^^1^2 * 0. As mentioned before, the solution 
becomes purely analytical for such cases. Figures 8 and 9 show the impeller forces 
at partial capacities with impeller and diffuser losses considered and neglected, 
respectively. For c o/Q < 0, F n increases and F t decreases with the decrease in 
$. This tendency can be understood from the discussions for a ♦ 0 in the last 
section through the effect of $ on v r 2 - The destabilizing region is enlarged by 
the decrease in <|>. It is interesting to note here that the value of w/Q 
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for which Imag(v r 2 ) = 0 increases with the decrease in <j>. In Figure 7, the 
impeller forces for 4>/4>d = 0*2 change largely near co/Q = 0.35. This condition 

corresponds to the rotating stall onset flow rate and the stall propagation velocity 
which can be determined by putting the denominator of equation (12) to be zero. It 
can be shown that the onset flow rate thus determined equals the flow rate at which 
the difference between the diffuser outlet static pressure and the impeller inlet 

total pressure becomes unchanged by a change in flow rate. Those conditions agree 

with those obtained by Railly & Ekerol (ref. 11). The onset condition determined 

from equation (12) for this case is (4>/4>d » = (0.21, 0.34). If we neglect the 

losses, as shown in Figure 8, we do not have the large change in the forces. For 
go/Q < 0, the losses have little effect on F n but they have the effect to increase 
F t near the design capacity and to decrease F t at lower capacities. This is 
because the change in losses due to the change in incidence angle becomes 

significant at lower capacities. In Figures 8 and 9, we observe destabilizing 
regions in co/Q < 0 at very low capacities. This has not been observed 

experimentally, and it is questionable if the present model can be applied to such a 
small capacity. 


Impeller forces at partial capacity 
(Cases with small impeller/diffuser gap) 

At small capacities the impeller forces change largely with a small change in 
r 3 / r 2 anc * the effect of rotating stall observed in Figure 7 for ** 3/^2 = 1 cannot 
be seen for ^ 3 /^ = 1.03. Figure 10 shows the impeller forces for r^/^2 = 1*10- 
The impeller forces change greatly in the region co/fi < 0 at partial capacities. 
The denominator of equation (12) becomes zero for (<t>/4>d » = (0.288, -0.16). 
Unlike the cases with ** 3 /^ = 1.0 shown in Figures 8 and 9, the large change does 
not disappear by neglecting impeller and diffuser losses. 

The vaned diffuser rotating stall onset condition can be obtained by a method 
similar to the present study under the assumptions that the pressure is constant at 
the outlet of the diffuser and that the flow upstream is irrotational . If we 
neglect the delay of loss, the onset flow rate coincides with that at which the 
difference of the inlet total pressure and the outlet static pressure is unaltered 
by the change in the flow rate under a constant inlet circumferential velocity. For 
the present case the diffuser rotating stall onset condition thus determined is 
(^/^d* w/Q) = (0.320, 0.0396). The impeller rotating stall onset condition for the 
present case determined under the condition of constant impeller outlet pressure is 
(4>/<t>d* w/Q) = (0.349, 0.914) (ref. 10). Both of these conditions are far apart from 
the condition above mentioned. 

As the flow rate gets smaller, the flow angle in the vaneless region becomes 
smaller and the time required for a fluid particle to pass through the region 
becomes larger. We consider a case when the shed vorticity decays proportionally to 
exp(-D a t), where Da is a positive constant and t is the time passed after the 
vorticity is shed from the impeller. Such cases can be treated by simply 
multiplying -D a ir(r^ - r 2 ^)/Q by \(r,0) in equation (4). Figure 11 shows the 
results with D a /Q = 2.0. We do not have the large change in co/Q < 0 as seen in 
Fig. 10. This fact may show that the large change is caused by the interaction of 
the vorticity with the vaned diffuser. The denominator of equation (4) becomes zero 
for (<t>/4>d» to/Q ) as 0.40, 0.02), and the impeller forces change greatly near this 
condition. This condition is close to that of diffuser rotating stall onset. This 
large change in the impeller forces disappears if we neglect the diffuser loss. 
Hence, the large change for D a /Q = 2.0 is considered to be related to the diffuser 
rotating stall. 
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Impeller forces at partial capacity 
(Cases with large impeller/diffuser gap) 


Figure 11 shows the results for 1 : 3 /^ = 1.5 and D a /Q = 0.25. The denominator 
of equation ( 12 ) becomes zero at ((J)/^, w/Q) = (0.33, 0 . 02 ) and (0.08, 0.98), near 
which the impeller forces change irregularly. The former and latter solutions 
disappear if we neglect the diffuser and impeller losses respectively. The diffuser 
and the impeller rotating stall onset conditions when they are alone are ( 4 >/<t>d» w/Q) 
= (0.328, 0.021) and (0.349, 0.914), respectively, and these are close to the above 
conditions. Hence the irregular changes of the forces in Figure 12 are considered 
to be related to the diffuser and impeller rotating stalls. The change of the 
forces near the diffuser stall onset condition is smaller than that in Figure 11 
for ^ 3/^2 = 1.1. It decreases further by the increase in r 3 /r 2 - On the contrary, 
the impeller stall onset flow rate increases with the increase in r 3 /r 2 , and its 
effect on the impeller forces becomes more significant. 

In any case the impeller forces change greatly near the diffuser and impeller 
rotating stall onset conditions and may become destabilizing for whirl. However, we 
should use a non-linear analysis for the exact determination of the impeller forces 
near rotating stall onset conditions, since rotating stalls are highly non-linear 
phenomena. In the present calculations with vaned diffusers, no phenomenon was 
observed which could be related to rotating stalls in the vaneless part of the 
diffuser . 


CONCLUSIONS 

Major findings of the present study may be summarized as follows. 

(1) The fluid forces on an impeller whirling in a vaned diffuser may become 
destabilizing for whirling motion even at design capacity. Discussions are made to 
elucidate the cause of the destabilizing fluid forces, and it was shown that the 
destabilizing fluid forces are closely related to the steady head-capacity 
characteristics of the impeller. 

(2) The destabilizing whirl velocity region is larger for a smaller diffuser vane 
angle, impeller/diffuser gap, and flow rate. 

(3) For an infinitesimally small gap between impeller and diffuser, the effects of a 
rotating stall on the impeller forces are observed at partial capacity. These 
effects disappear with a small increase in the gap, but the effects of the 
interaction of shed vorticity with the vaned diffuser appear instead. If we take 
into account the decay of shed vorticity, effects of diffuser rotating stalls are 
observed. 

(4) As we increase the gap further, the effects of the diffuser rotating stall 
become smaller. The impeller rotating stall onset flow rate increases and its 
effect on the impeller forces becomes more significant. In either case the impeller 
forces change greatly near the rotating stall onset conditions and may become 
destabilizing for impeller whirl. 

(5) For small positive to negative whirl velocity ratio, the contribution of the 
impeller outlet pressure distribution has the most significant effect on the impeller 
forces. It increases as the impeller/diffuser gap is decreased. 
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APPENDIX 

(Normal and tangential forces for ^^ 1^2 * 1 and ot + 0, without loss) 


From equation (14) and (15) we obtain 
F t = (M/r2)(l - u/Q)<t> 

F n = (M/r2)(l - u/fi)('|i - w/Q) 

where <|> r Q/ 2 irr 2 ^ Q is the flow coefficient and t|» = 1 - <{> cot (3 is the head 
coefficient (normalized by p^^Q^) of the impeller. M is the inertial length of 
the impeller flow channel, defined in the paragraph following equation (6). 

These expressions show that Fj- is linear with respect to co/Q and is 
proportional to the flow coefficient ♦ • F n is a parabolic function of u/Q and 
dependent on the head coefficient These relations are illustrated in the 

following figure. 



Illustration of Ft and F n for the cases with rz/vz ■+ 1 and a *+ 0, without losses. 
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Fig. 4 Effects of diffuser vane angle on whirling 
forces at design capacity. The effective 
length of diffuser flow channel is kept 
constant. 



Fig. 5 Division of whirling forces. Contribution of 
impeller outlet pressure distribution is 
dominant in a range of small positive to nega- 
tive whirl velocity ratio. 



Fig. 6 Illustration of the relations among impeller 
outlet flow distribution, pressure distribu- 
tion and impeller forces, for an idealized 
case with a + 0 and ra •* ra. 



Fig. 7 Illustration of the relations among impeller 
outlet flow distribution, pressure distribu- 
tion and impeller forces, for a realistic 
case. 
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With imp. 8, diffu. losses 


r 3 /r 2 =1.0 
r t /rj=1.25 
a = 1 1 deg. 
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Fig. 8 Impeller forces at partial capacities, for 
the case with infinitesimally small impeller/ 
diffuser gap, with the effects of impeller and 
diffuser losses taken into account. Effects 
of a rotating stall are observed. 


Fig. 9 Impeller forces at partial capacities, for 

the case with infinitesimally small impeller/ 
diffuser gap. The impeller and diffuser 
losses are not considered. We do not observe 
the effects of rotating stall. 


With imp. & diffu. losses 


With imp. & diffu. losses 


r L l r 2 =1.35 
a sll deg. 
Da/fl =0.0 


r 3 /r 2 = 1.10 
rj r 2 = 1.35 
a = 1 1 deg 
Da/O=2.0 


Fig. 10 Impeller forces at partial capacities for 
the case with smaller impeller/diffuser gap. 
We observe fluctuations of the forces in the 
range of negative whirl velocity. 


Fig. 11 Impeller forces at partial capacities for 
the case with smaller impeller/diffuser gap. 
Exponential decay of shed vorticity is 
considered. Fluctuations in the range of 
negative whirl velocity diminish and the ef- 
fects of diffuser rotating stall are observed at 
small positive whirl velocity. 












